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j  '  Improved  spectral  estimation  techniques  are  becoming  a  valuable  asset  in  adaptive  processing  array 
antenna  systems.  They  are  invaluable  because  of  the  amount  of  additional  useful  information  they  can 
provide  about  the  interference  environment,  utilizing  only  a  relatively  small  number  of  degrees-of- 
freedom  (DOF).  The  •superresolution'’  capabilities,  estimation  or  coherence,  and  relative  power  level 
determination  serve  to  complement  and  refine  the  data  from  faster  conventional  estimation  techniques.  •  v. 

This  report  discusses  two  conceptual  application  area  examples  that  use  such  techniques;  partially- 
adaptive  low-sidelobe  arrays,  and  fully-adaptive  tracking  arrays.  For  the  partiaily-adaptive  area,  the  infor¬ 
mation  is  used  for  the  efficient  assignment  of  a  limited  number  of  DOF  in  a  beamspace  constrained  adap¬ 
tive  system  to  obtain  a  stable  mainbeam,  retention  or  low  sidelobes,  considerably  faster  response,  and 
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19.  ABSTRACT  (Continued) 

reduction  in  overall  cost.  These  benetits  arc  demonstrated  via  simulation  examples  computed  for  a  16- 
elenient  linear  array.  For  the  fully  adaptive  tracking  airay  area,  the  information  is  used  in  an  all-digital 
processing  system  concept  to  permit  stable  nulling  of  coherent  interference  sources  in  the  mainbeam 
region,  efficient  assignmeni/conirol  of  the  available  DOF,  and  greater  flexibility  in  time-domain  adaptive 
filtering  strategy. 


CONTENTS 


INTRODUCTION  . -  1 

PARTiALLY-ADAPTIVE  LOW-S1DELOBE  ANTENNAS  .  1 

Low-Sidelobe  Eigenvector  Constraint  .  2 

Low-Sidelobe  Constraints  tor  a  General  Beamlormer  .  6 

Interference  Sources  in  the  Mainbeam  Region  .  9 

SOURCE  ESTIMATION  AND  BEAM  ASSIGNMENT  .  9 

AN  ADAPTIVE  ARRAY  TRACKING  APPLICATION  .  16 

Coherent  Spatial  Interference  Sources  .  16 

All-Digital  Tracking  System  Concept  .  18 

CONCLUSIONS .  20 

REFERENCES .  21 

APPENDIX  A  —  Snapshot  Signal  Model  .  23 

APPENDIX  B  —  Eigenvalue/ Eigenvector  Decomposition .  26 

APPENDIX  C  —  Tracking  Beams  and  Adaptive  Distortion  .  33 

APPENDIX  D  —  Forward-Backward  Subaperture  Averaging  .  37 


iii 


USING  SPECTRAL  ESTIMATION  TECHNIQUES  IN  ADAPTIVE 
PROCESSING  ANTENNA  SYSTEMS 

INTRODUCTION 

Improved  spectral  estimation  techniques  are  an  emerging  technology  which  derives  largely  from 
modern  spectral  estimation  theory  of  the  past  decade  and  adaptive  array  processing  techniques  [1,2,3] 
Coupled  with  the  phenomenal  advances  in  digital  processing,  these  techniques  are  becoming  a  valuable 
asset  for  adaptive  array  antenna  systems.  Their  value  lies  in  the  considerable  amount  of  additional  use¬ 
ful  information  which  they  can  provide  about  the  environment  while  using  only  a  relatively  small 
number  of  degrees-of-freedom  (DOF).  For  example,  current  spectral  estimation  algorithms  can  pro¬ 
vide  asymptotically  unbiased  estimates  of  the  number  of  interference  sources,  source  directions,  source 
strength  ,  and  any  cross-correlations  (coherence)  between  sources  [4,5].  Such  information  can  then  be 
used  to  track  and  "catalogue"  interference  sources,  hence  assign  adaptive  DOF. 

These  newer  techniques  are  not  viewed  as  a  "superresolution"  replacement  for  more  conventional 
estimation  methods  such  as  mainbeam  search,  analogue  beamformers,  or  spatial  discreet  Fourier 
transforms  (DFT);  but  rather,  the  new  technology  is  considered  complementary  to  the  other  methods 
and  best  used  in  tandem.  For  example,  "superresolution"  techniques  cannot  compete  with  the  speed  of 
a  DFT  Some  comparisons  of  the  various  methods  may  be  found  described  in  the  literature  (3,5,6). 

The  purpose  of  this  report  is  to  present  two  conceptual  application  areas  that  use  spectral  estima¬ 
tion  techniques,  partially-adaptive  low-sidelobe  antennas,  and  fully-adaptive  tracking  arrays.  In  a  par¬ 
tially  adaptive  array  only  a  part  of  the  DOF,  array  elements  or  beams,  is  individually  controlled  adap¬ 
tively  [7,8,9].  Obviously,  the  fully  adaptive  configuration  is  preferred  since  it  offers  the  most  control 
over  the  response  of  the  antenna  system.  But,  when  the  number  of  elements  or  beams  becomes 
moderately  large  (hundreds),  the  fully  adaptive  processor  implementation  can  become  prohibitive  in 
cost,  size,  and  weight. 

This  report  is  divided  into  three  principal  parts.  The  second  section  of  this  report  discusses 
partially-adaptive,  low-sidelobe  antennas  with  the  focus  upon  a  constrained  beamspace  system.  In  the 
third  section  of  the  report,  the  source  estimation  and  beam  assignment  from  superresolution  techniques 
are  considered;  and  in  the  fourth  section,  an  all-digital,  fully-adaptive  tracking  array  concept  is  dis¬ 
cussed.  Several  Appendices  are  also  used  for  referral  of  the  more  tedious  details. 

PARTIALLY-ADAPTIVE  LOW-SIDELOBE  ANTENNAS 

The  antenna  system,  addressed  in  this  section,  is  assumed  to  be  a  moderately  large  aperture  array 
of  low-sidelobe  design  wherein  the  investment  is  already  considerable  and  one  simply  could  not  afford 
to  make  it  fully  adaptive.  The  assumption  of  low-sidelobes  (30  dB  or  better)  is  intended  to  give  us 
good  initial  protection  against  modest  interference  sources  and  to  reduce  the  problems  from  strong 
sources,  i.e.,  in  regard  to  the  number  of  adaptive  DOF  required  and  the  adaptive  dynamic  range  of  the 
processor.  Thus,  retention  of  the  low  sidelobes  is  considered  a  major  goal  in  our  adaptive  system.  In 
the  discussion  that  follows  it  is  shown  that  using  improved  spectral  estimation  techniques  in  such  a  sys¬ 
tem  can  result  in  the  following  benefits  over  a  fully  adapthe  array  system: 
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a.  Reduction  in  overall  cost  because  relatively  few  adaptive  DOF  are  implemented. 

b.  Simple  adaptive  weight  constraints  permit  minimal  degradation  of  the  mainbeam  and  sidelobe 
levels. 

c.  Reduction  in  computation  burden. 

d.  Considerably  faster  adaptive  response. 

e.  Compatibility  with  a  larger  number  of  adaptive  algorithms,  including  analogue  versions. 

f.  Greater  flexibility  in  achieving  a  "tailored1'  response  due  to  greater  information  available. 

On  the  negative  side,  a  partially-adaptive  system  can  never  be  guaranteed  a  cancellation  perfor¬ 
mance  equal  to  that  of  a  fully  adaptive  array,  and  will  deteriorate  abruptly  in  performance  when  the 
interference  situation  exceeds  its  adaptive  DOF.  These  risks  are  an  inherent  pait  of  the  package  and 
must  be  carefully  weighed  for  any  specific  system  application. 


Low-Sidelobe  Eigenvector  Constraint 


In  this  section  we  review  how  unconstrained  adaptive  arrays  can  experience  very  "noisy"  sidelobe 
fluctuations  and  mainbeam  perturbations  when  the  data  observation/integration  time  is  not  iong 
enough,  even  though  the  quiescent  mainbeam  weights  are  chosen  for  low  sidelobes.  Consider  the  sim¬ 
ple  schematic  shown  in  Fig.  1,  and  let  us  compute  the  complex  adaptive  element  weights  Wk  from  the 
well-known  Sample  Matrix  Inverse  (SMI)  algorithm  [9,10]  as  expressed  in  the  following  matrix  nota¬ 
tion. 


where: 


w~mr-'s* 

W  is  the  adaptive  weights  vector, 

R  is  the  sample  covariance  matrix, 

S*  is  the  quiescent  mainbeam  weights  vector,  and  n  is  a 
constant. 


(1) 


*  denotes  the  conjugate  of  a  complex  vector  or  matrix. 


F,g  1  —  Schematic  of  adaptive  urray 


Furthermore,  the  sample  covariance  matrix  is  computed  via  the  simple  "block"  average  taken  over  N 
snapshots. 


(2) 
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where  E («)  is  the  element  signal  data  vector  received  at  the  nth  time  sampling.  (See  Appendix  A  for 
the  description  of  snapshot  signal  model.)  The  data  observation/integration  time  in  (2)  is  the  parame¬ 
ter  N.  If  R  is  estimated  over  a  lengthy  observation  time,  like  thousands  of  snapsh  >ts,  then  the  sidelobe 
fluctuations  from  W  updates  will  be  relatively  small.  However,  practical  system  usage  often  demands 
short  observation  times  on  the  order  of  hundreds  of  snapshots  or  even  less. 

Figure  2  illustrates  a  typical  adapted  pattern  behavioi  for  independent  estimates  of  R  using  N  — 
256  snapshots  per  update  for  the  case  of  three  30  dB  noncoherent  sources  located  at  14,  18  and  22°. 
The  antenna  aperture  chosen  for  this  example  is  a  16-element  linear  array  with  half-wavelength  element 
spacing  and  a  30  dB  Taylor  illumination  incorporated  in  S’.  The  adaptive  algorithm  maintains  the 
mainbeam  region  and  successfully  nulls  out  the  interference  sources,  bui  it  also  raises  the  sidelobe  lev¬ 
els  elsewhere.  The  adaptive  patterns  are  in  continual  fluctuation  in  the  sidelobe  regions  and  may 
exceed  the  quiescent  sidelobe  level  by  a  considerable  margin.  Also,  the  mainbeam  suffers  a  significant 
modulation  which  would  degrade  tracking  performance.  These  effects  worsen  as  the  value  of  N 
decreases. 

To  understand  the  reason  for  this  undulating  pattern  behavior,  it  is  helpful  to  analyze  the 
optimum  weights  in  terms  of  eigenvalue/eigenvector  decomposition.  Appendix  B  contains  a  derivation 
of  such  a  decomposition  for  Eq.  (1),  and  we  reproduce  Eq.  CB- 18)  below: 


W 


where: 


K 

S’-  I 

a!  -  Po 

a,e, 

a7 

'-1 

p i 

-  _  ..  */f  * 


(3) 


i  denotes  the  transpose  of  a  vector  or  matrix.  The  /3,2  and  e,  are  the  eigenvalues  and  eigenvectors, 
respectively,  of  the  sample  covariance  matrix,  and  ,flo  is  equal  to  receiver  channel  noise  power  level. 
Equation  (3)  shows  that  VV  consists  of  two  parts:  the  first  pari  is  the  quicsceni  mainbeam  weight  S’; 
the  second  part,  which  is  subtracted  from  S’,  is  the  summation  of  weighted  orthogonal  eigenvectors. 
This  is  a  clear  expression  of  the  fundamental  principle  of  pattern  subtraction  which  applies  in  adaptive 
array  analysis  19,1 1]. 


We  introduce  the  term  principal  eigenvectors  (PE)  to  mean  those  eigenvectors  which  correspond  to 
unique  eigenvalues  generated  by  the  spatial  source  distribution;  and  the  term  noise  eigenvectors  to  mean 
those  eigenvectors  which  correspond  to  the  small  noise  eigenvalues  generated  by  the  receiver  channel 
noise  contained  in  the  finite  R  estimates.  The  PE  are  generally  rather  robust  arid  tend  to  remain  rela¬ 
tively  stable  from  one  data  trial  to  the  next,  whereas  the  noise  eigenve*  ors  tend  to  fluctuate  consider¬ 
ably  because  of  the  inherent  random  behavior  of  noise.  This  difference  in  behavior  is  illustrated  in  Fig. 
3  for  the  three  source  case  described  above;  wherein  there  are  three  PE  and  thirteen  noise  eigen*  actors 
associated  with  each  R  estimate.  Figure  3(a)  shows  the  stability  of  the  three  PE  for  nine  trials,  and  Fig. 
3(b)  shows  the  random  behavior  of  typical  noise  eigenvectors  for  the  exact  same  trials.  Thus,  we 
would  expect  that  the  sidelobe  undulations  in  Fig.  2(b)  are  associated  primarily  with  the  noise  eigenvec¬ 
tors.  This  thesis  is  verified  in  Fig.  4,  which  illustrates  the  adapted  patterns  resulting  from  Eq.  (3)  when 
only  the  PE  are  subtracted 


The  above  adaptive  array  pattern  behavior  leads  to  the  following  observations  for  source  distribu¬ 
tions  which  do  not  encroach  upon  the  mainbeam  and  involve  a  small  number  of  the  available  degrees- 
of-freedom  (DOF): 

a.  It  is  possible  to  retain  low  sidelobes  in  the  adapted  patterns,  even  with  short  observation  times, 
by  constraining  our  algorithm,  Eq.  (3),  to  utilize  oniy  the  PE.  The  weight  solution  is  unique  and  there¬ 
fore  stable. 
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SPATIAL  ANGLE  IN  DEGREES 

Fig  4  —  Typical  adapted  patterns  resulting  from  the  constraint 
of  utilizing  only  the  PL,  three-source  case  of  Fig  2 


b.  Utilizing  only  the  PE  is  tantamount  to  operating  our  adaptive  system  in  beamspace  (as  opposed 
to  element  space)  with  a  set  of  weighted  orthogonal  canceller  beams. 

c.  The  fully  adaptive  array  automatically  forms  and  "assigns”  its  PE  canceller  beams  to  cover  the 
interference  source  distribution,  with  one  beam  per  each  DOF  needed. 

Therefore,  we  have  set  forth  a  low-sitlelobe  eigenvector  constraint  algorithm  for  this  type  of  res¬ 
tricted  interference  situation. 


Low- Sideiobe  Constraints  for  a  General  Beaniformer 

Consider  next  a  more  interesting  configuration  which  is  shown  by  the  schematic  diagram  in  Fig. 
5,  where  we  represent  an  adaptive  array  system  operating  in  beamspace  so  as  to  have  available  some 
pre-adaptation  spatial  filtering.  Applebaum  and  Chapman  [8,9,12]  were  the  first  to  describe  beamspace 
systems  of  this  type,  using  a  Butler  matrix  1 1 3]  beamformer,  wherein  the  vector  of  the  beamformer 
outputs,  E,  may  be  expressed  as  follows: 

E-B'E  (4) 

where  B  is  a  KxK  matrix  containing  the  beamformer  element  weights  (see  Appendix  C).  Other 
descriptions  of  beamspace  systems  are  also  available  in  the  literature  [9,14,15,16],  of  which  Adams  et 
al.  [15],  is  particularly  germane  to  our  discussion.  Chapman  [8]  pointed  out  that  when  used  in  a 
partially  adaptive  configuration,  such  beamspace  systems  are  susceptible  to  aperture  element  errors  and 
cannot  arbitrarily  compensate  the  random  error  component  of  their  sideiobe  structure.  This  makes  it 
necessary  to  control  element  errors  in  accordance  with  the  quiescent  mainbeam  sideiobe  level  desired, 
and  fits  into  our  initial  assumption  of  low-sidelobe  design  as  mentioned  earlier.  A  separate  weighted 


6 


•?o 


NRL  RLP0H1  8920 


mainbeam  summing  is  indicated  which  may  be  obtained  either  by  coupling  into  the  beamformer  out¬ 
puts  as  shown,  or  by  coupling  off  from  the  elements  and  providing  suitable  phase  shifters  for  steering 
plus  a  corporate  feed  network.  Our  purpose  here  is  to  examine  the  sidelobc  performance  of  such  a 
partially-adaptive  beamspacc  system  in  which  element  errors  are  kept  low  and  beamformci  beams  are 
subjected  to  simple  constraints. 

Spatiai  estimation  data  on  the  interference  source  distribution  shall  determine  which  beamformer 
beams  are  to  be  adaptively  controlled.  Such  beams  are  defined  herein  as  assigned  beams,  and  the  idea 
is  to  assign  only  enough  beams  to  accommodate  the  DOF  required  by  the  source  distribution.  When¬ 
ever  the  two  are  equal,  the  adaptive  weight  solution  is  unique  and  we  avoid  adding  any  extra  "noisy" 
weight  perturbations.  The  reader  will  recognize  ihat  we  are  attempting  to  replace  the  PE  beams  of  the 
previous  section,  Lcw-Sidelobe  Eigenvector  Constraint,  with  assigned  beams  from  our  general  beam- 
former.  Thus,  we  are  defining  a  partially-adaptive  array  which  will  utilize  only  a  lelatively  small 
number  of  its  available  DOF.  In  addition  to  this  assigned  beam  constraint,  we  seek  to  limit  the  adap 
tive  weights  of  assigned  beams  to  a  maximum  level,  y,  which  was  chosen  to  exceed  the  mainbeam 
sidelobe  level  by  only  a  few  dB.  This  prevents  an  excessive  rise  ir  adaptive  sidelobe  level,  including 
the  condition  where  the  number  of  assigned  beams  exceeds  the  DOE  required. 

y  represents  the  product  of  assigned  beam  gain  and  adaptive  weight  magnitude,  so  that  we  have 
the  option  of  working  with  beamformer  beams  which  are  considerably  decoupleii/attenuated. 

An  equation  formulation  may  be  expressed  in  terms  of  the  same  pattern  subtraction  principle  as 
used  in  Eq.  (3)  for  K  beams, 

W„-S*-  £  (5) 

k~\ 
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whcrc: 

I  Wk  I  <  y  lor  J  assigned  beams 
—  0  for  all  other  beams 

and  b*  is  the  A th  Butler  matrix  beam  clement-weight  vector.  When  Wk  ~  0,  that  beam  port  is  essen¬ 
tially  disconnected  from  the  output  summation  and  it  is  much  to  our  advantage  to  reduce  the  DOF  of 
the  adaptive  weight  processor  accordingly,  i.e.,  this  processor  reduction  relates  directly  to  the  computa¬ 
tional  burden,  response  time,  sideiobc  degradation,  and  overall  cost  mentioned  earlier.  For  example, 
utilizing  the  SMI  technique  described  in  Eqs.  (1)  and  (2),  we  would  now  hav*'  the  advantage  that  our 
sample  covariance  matrix  of  signal  inputs,  R,  involves  only  the  J  assigned  beams  and  its  dimensions 
reduced  from  KxK  to  JxJ,  thereby  greatly  easing  the  computation  burden  involved  in  obtaining  its 
inverse  (9).  The  equivalent  "steering  vector",  A,  per  Applebaum  [7]  is  also  reduced  to  dimension  J  and 
consists  of  the  cross-correlation  between  the  mainbeam  signal  V  and  the  J  assigned  beam  outputs,  Y, 

A  -  4:  y  F  ( u  ‘  Y  * ( n )  .  (6) 

N  4*i 

The  y th  assigned  beam  output  for  the  nth  snapshot  signal  sample  is  simply 

Yj(n)  —  E'(n)bA  ,  k  set  by  j  (7) 

where  the  particular  value  of  k  must  be  selected  for  the  ylh  assigned  beam  Our  J  dimension  adaptive 
weight  solution  thus  becomes, 

W  -  R  'A.  (8) 

Equation  (8)  gives  us  the  J  assigned  beam  weights  required  in  Eq.  (5).  The  proposed  constraint 
I  Wk  I  <  y  can  be  applied  directly  to  the  solution  from  (8),  but  be  aware  that  this  is  a  "hard"  constraint 
and  the  results  will  not  be  optimal  when  the  limit  is  exceeded. 

A  softer,  more  flexible  constraint  for  our  purposes  is  one  suggested  by  Brennan*  based  upon  Ows¬ 
ley  [17],  where  weights  are  selected  which  simultaneously  minimize  botii  the  output  and  the  sum  of  the 
weight  amplitudes  squared,  i.e., 

minimize  |  V  -  W'Yl2  +  aW'W* 

where  the  overbar  denotes  averaging  over  N  snaps.  The  solution  is  a  simple  modification  to  Ea.  (8) 
wherein 

W  -  |  R  +  al  ]  '  A  (9) 

where: 

a  -  -*j-  Trace  [R] 

We  note  here  that  Eq.  (9)  adds  a  small  percentage  of  the  average  assigned  beam  power  to  the  diagonal 
terms  of  R  We  also  recall  that  >■  was  selected  to  be  close  to  the  mainbeam  sideiobc  level.  Although  a 
is  a  small  percentage  of  the  Trace  [R],  it  is  generally  much  larger  than  the  receiver  noise  level,  /J(?;  this 
domination  over  receiver  noise  by  a  constant  will  tend  to  sevciely  dampen  weight  fluctuations  due  to 
noise.  Equation  (9)  deviates  from  the  optimum  Weiner  weights  and  will  result  in  a  slightly  larger  out¬ 
put  residue,  however,  the  cost  is  negligible  compared  to  the  remarkably  stable  resuits  achieved  from 
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this  rather  simple  constraint.  It  essentially  permits  the  number  of  assigned  beams  to  exceed  the  DOF 
required,  and  yet  retain  low  sidelobe  levels. 

Equations  (S)  through  (9)  were  used  in  computing  the  adaptive  pattern  examples  which  follow.  The 
reader  should  recognize  that  the  J  dimension  adaptive  weight  solution  may  be  arrived  at  via  any  of  the 
current  adaptive  processing  algorithms  such  as  Howells- Applebaum  [7] ,  Gram-Schmidt  [9),  Sample 
Matrix  Inverse  Update  (18],  etc. 

When  applying  these  constraints  to  our  three-source  case  of  Fig.  2,  we  would  assign  beamformer 
beams  Nos.  10,  11,  and  12  to  cover  the  sources  as  illustrated  in  Fig.  6(a).  These  assigned  beams  are 
then  given  a  maximum  gain  level  of  about  5  dB  above  the  -30  dB  mainbeam  sidelobes;  thus,  the 
assigned  beam  weights  are  constrained  to  I  Wk  |  <  0.055.  All  other  Wk  are  set  to  zero.  Typical  resul¬ 
tant  adapted  patterns  are  shown  in  Fig.  6(b),  where  nine  trials  of  160  snapshots  each  are  plotted.  The 
pattern  stability  is  near-perfect  for  a  unique  solution  like  this,  and  note  that  the  three  sources  have 
been  nulled  with  very  little  perturbation  of  the  mainbeam  sidelobes  except  in  the  immediate  vicinity  of 
the  sources.  Since  we  are  inverting  a  matrix  of  only  3x3  dimension  in  Eq.  (8),  for  this  case,  it  follows 
that  the  number  of  snapshots  processed  per  trial  could  be  reduced  by  an  order  of  magnitude  [10]  and 
still  obtain  excellent  results. 

Figure  7  demonstrates  how  the  adaptive  weights  will  become  noisy  if  we  include  even  one  extra 
DOF  beyond  the  unique  solution.  Beamformer  beam  No.  16  was  deliberately  added  for  the  same  case 
as  in  Fig.  6,  and  we  may  note  the  consequent  sidelobe  fluctuations.  However,  if  we  use  the  "soft"  con¬ 
straint  of  Eq.  (9)  in  solving  for  the  weights,  stable  performance  is  again  restored  despite  the  extra  DOF. 
It  may  be  of  interest  to  the  reader  that  for  this  particular  example: 

J  -  4 

y  <*  0.045  (-27  dB) 

Trace  [R)/J  s=  10,350  /302 
«  -  21/3o  • 

Although  not  shown  here,  another  example  of  interest  is  the  case  of  using  a  two-beam  cluster 
(Nos  1 1  and  12)  to  cancel  a  single  40  dB  broadband  source  located  at  22°.  It  was  found  that  the  source 
could  be  adequately  cancelled  at  bandwidths  up  to  15%. 

Many  other  .•->mbinations  of  source  distributions  and  assigned  beams  were  tested  to  further  verify 
the  teennique,  arid  the  partially-adaplive  performance  was  satisfactory  provided  that  the  assigned  beams 
were  sufficient  to  cover  the  DOF  demanded  by  the  source  distribution 

Interference  Sources  in  the  Mainbeam  Region 

Extension  of  the  foregoing  partially-adaptive  array  technique  for  mainbeam  interference  is 
straightforward,  provided  we  relax  the  constraint  upon  the  value  of  y  in  Eq.  (5).  Obviously,  the  low- 
sidelobe  stratagem  becomes  secondary  to  the  greater  menace  of  an  interference  source  coming  in 
through  our  high-gain  mainbeam.  Low  sidelobes  could  still  be  retained,  if  necessary,  by  implementing 
a  beamformer  which  is  capable  of  producing  a  family  of  low-sidelobe  assigned  beams  [15]. 

SOURCE  ESTIMATION  AND  BEAM  ASSIGNMENT 

Modern  spectral  estimation  techniques  are  a  welcome  addition  to  the  conventional  methods  for 
tracking  and  cataloging  interference  sources.  They  do  not  interfere  with  any  functions  of  the  main- 
beam,  and  are  capable  of  providing  superior  source  resolution  from  fewer  elements.  The  latter  advan¬ 
tage  is  obtained  in  part  because  we  assumed  low  sidelobes  for  the  mainbeam,  i.e.,  the  only  sources  that 
require  estimation  are  those  few  which  are  of  sufficiently  high  SNR  to  get  through  the  mainbeam 
sidelobes.  Resolution  performance  is  always  directly  related  to  SNR,  [2,5,6], 
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The  principle  ror  achieving  source  estimation  from  a  small  fraction  of  the  aperture  X)F  has  been 
demonstrated  via  many  techniques,  both  conventional  and  optimal  (1,3,19).  It  is  not  wuhin  the  scope 
of  this  report  to  attempt  a  comprehensive  comparison  of  such  techniques,  but  the  point  is  important  to 
our  Concept  so  that  an  example  of  a  half-aperture  linear  array  estimator  is  given  in  this  section.  The 
type  of  application  envisioned  is  illustrated  in  Fig.  8,  where  we  represent  a  KxK  element  aperture  sys¬ 
tem  in  which  the  adaptive  beam  DOF  are  to  be  assigned  on  the  basis  of  estimates  derived  from  two 
orthogonal  linear  arrays  of  Kj 2  elements  each.  An  extension  of  the  2D  (two-dimension)  beamspace 
adaptive  array  system  of  Fig.  5  to  the  3D  system  suggested  by  Fig.  8  permits  several  beamformer 
options,  including: 
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Fig.  8  —  (K  x  K)  element  apertuie  within 
which  row/column  linear  arrays  couple  into 
source  estimation  processors 


Fig.  9  —  Typical  RF  receiver  techniques  associated 
with  A/D  complex  data  sampling 
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a.  Two  orthogonal  2D  beamformcrs  of  which  one  is  coupled  into  a  row  and  the  other  coupled 
into  a  column  of  elements. 

b.  A  complete  3D  beamformer  [201  coupled  into  the  aperture  elements,  perhaps  on  a  thinned 

basis. 

The  separate  mainbeam  must  be  summed  from  all  A'2  elements  to  attain  the  desired  low  sidelobes. 

Although  they  involve  relatively  few  elements  from  the  aperture,  the  linear  array  estimators 
represent  a  significant  increase  in  system  expense,  because  they  are  all-digital  processing  subsystems. 
Typical  RF  receiver  components  required  prior  to  the  signal  analogne-to-digital  (A/D)  converters  are 
shown  in  Fig  9.  The  processing  of  the  digital  signals  to  estimate  the  sources  may  be  carried  out  in 
accordance  with  a  number  of  spectral  estimation  algorithms  available  in  the  literature  [1-6].  Appendix 
B  discusses  several  algorithms  that  were  used  in  the  simulations  conducted  for  this  report.  For  exam¬ 
ple,  Fig.  10  illustrates  a  comparison  plot  of  our  mainbeam  search  scan  vs  half-aperture  eigenanalysis 
processing  results  for  the  16-element  linear  array  case  of  Fig.  2. 


spatial  angle  in  degrees 


Fig.  10  —  Comparison  of  mainbeam  scan  vs  half-aperture  eigenanalysis 
source  estimation  for  three  source  case  of  Fig.  2 

Once  the  source  estimation  information  is  available,  we  can  proceed  by  assigning  beamfoimer 
beams  via  a  computer  logic  program.  For  the  simulations  reported  in  this  report,  a  Fortran  IV  com¬ 
puter  code  named  "BEAMASSIGN"  was  developed  which  accepts  source  information  updates,  compares 
the  new  data  against  a  source  directory  kept  in  memory,  computes  track  updates  for  sources  already  in 
memory,  determines  priority  ranking,  and  assigns  beams  to  cover  the  sources  of  nighest  priority.  An 
important  point  to  note  is  that  beam  assignment  does  not  require  great  accuracy,  i.e.,  a  half-beamwidth 
is  usually  close  enough.  Also,  clusters  of  two  or  three  adjacent  beams  may  be  assigned  for  doubtful 
cases. 
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A  demonstration  of  beam  assignment  was  conducted  with  a  moving  source  simulation  involving 
the  16-element  linear  array  of  Fig.  2.  Four  sources  of  unequal  strength  were  set  up  in  the  farfield, 
traveling  in  criss-crossing  patterns.  Two  of  the  sources  are  of  30  dB  strength  with  start  angles  of  3.0 
and  39. 0°,  and  two  are  of  43  dB  strength  with  start-angles  of  5.0  and  70.0°.  The  estimation  of  the 
scanned  mainbeam  for  this  example  is  shown  in  Fig.  11(a).  Each  time-unit  plot  cut  is  computed  from 
R  averaged  over  160  snapshots, 

PQ  =  S'RS‘  (10) 

where  S’  is  the  mainbeam  steering  vector  used  to  generate  the  display  plot.  As  expected,  this  simple 
Fourier  output  is  dominated  by  the  two  stronger  sources.  Ip.  contrast.  Fig.  11(b)  shows  the  source  esti¬ 
mation  derived  from  eigenanalysis  processing  using  only  half  of  the  aperture  (8  elements).  Note  that 
the  "superresolution"  characteristics  of  this  type  of  optimal  estimation  produces  excellent  source  track¬ 
ing,  even  in  the  vicinity  of  cross-over  of  three  of  the  sources. 


Fig  1 1  —  bstimation  of  four  moving  sources  via  mainbeam  scan  and  half  aperture  eigenanalysis  algorithm 
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The  results  from  the  use  of  the  source  information  data  contained  in  Fig.  11(b)  to  continuously 
update  beam  assignments  are  illustrated  in  the  adapted  pattern  cuts  shown  in  Fig.  12(a).  Note  that  the 
mainbc  tm  remains  steady,  and  the  sidelobcs  seldom  exceed  their  quiescent  30  dB  peak  level,  despite 
the  drastic  shifting  of  the  nulls  as  the  moving  sources  criss  cross  in  the  sidelobe  region  In  contrast, 
Fig.  12(b)  illustrates  the  adapted  pattern  cuts  obtained  when  we  utilize  the  SMI  algorithm  weights  with 
the  array  fully  adaptive.  Although  the  source  cancellation  is  excellent,  the  mainbeam  suffers  significant 
modulation  and  the  peak  sidelobe  levels  rise  considerably. 


Fig.  12  —  Adaptive  patterns  for  IS  element  linear  array.  SMI  algorithm,  four  moving  sources  case  of  Fig.  1 1 
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AN  ADAPTIVE  ARRAY  TRACKING  APPLICATION 

A  second  area  where  spectral  estimation  techniques  can  provide  valuable  assistance  is  that  of 
adaptive  array  tracking  systems.  Here  we  are  dealing  with  the  problem  of  attempting  to  track  targets 
under  the  condition  of  having  inteiference  sources  present  in  the  mainbcam  region.  Some  early  pro¬ 
posed  solutions  in  this  area  evolved  from  the  growing  adaptive  an  ay  technology  of  the  1970’s.  For 
example,  a  paper  by  White  (21)  discusses  the  radar  problem  of  tracking  targets  in  the  low-angle  regime, 
where  conventional  tracking  radars  encounter  much  difficulty  because  of  the  presence  of  a  strong 
surface-reflected  ray. 

The  first  extension  of  fully  adaptive  arrays  tc  angle  estimation  in  external  noise  fields  is  the  con¬ 
tribution  of  Davis  et  at.  (22),  who  developed  an  algorithm  based  on  the  outputs  of  adaptively  distorted 
sum  and  difference  beams.  The  adaptive  beams  filter  (null)  the  external  noise  sources,  and  distortion 
correction  is  then  applied  in  the  resultant  monopulse  output  angle  estimate.  Their  work  is  particularly 
appropriate  as  a  starting  point  for  this  section,  where  we  discuss  the  advantages  of  using  spectral  estima¬ 
tion  techniques  in  an  all-digital,  fully  adaptive,  array  tracking  system.  Reference  (15)  is  also  pertinent. 

Coherent  Spatial  Interference  Sources 

The  existence  of  significant  coherence  between  spatial  sources  as,  for  example,  in  multipath  situa¬ 
tions  involving  a  specular  reflection,  continues  to  represent  a  serious  problem  area  even  for  a  fully 
adaptive  tracking  array. 

Reasons  include: 

a.  Coherent  signals  in  space  are  not  stationary  (2,5,23). 

b.  Adaptive  systems  may  perform  cancellation  via  weight  phasing  rather  than  null  steering 
(5,23,24,25,26). 

c.  Adaptive  tracking  beam  distortion  is  highly  sensitive  to  coherent  signal  phasing. 

d.  Signal  fading  under  anti-phase  conditions. 

To  demonstrate  these  reasons,  adaptive  characteristics  were  computed  for  a  16-element  linear  "ray  for 
an  interference  case  in  which  there  are  two  13  dB  coherent  sources  in  the  mainbeam  region.  There  is 
also  a  third  source,  non-coherent,  in  the  nearby  sidelobe  region  that  acts  as  a  stable  null  comparison 
point. 

In  Fig.  13(a),  we  illustrate  the  severe  changes  in  our  mainbeam  caused  by  variation  of  the  phase 
shift  between  the  two  coherent  sources.  The  quiescent  mainbeam  has  the  same  Taylor  weighting  as 
that  in  Fig.  2(a).  Figure  13(b)  illustrates  the  spatial  insertion  loss  associated  with  the  three  adaptive 
weightings  invoived.  Note  that  for  source  phasing  of  0°  and  180°,  the.  adaptive  weights  do  not  achieve 
cancellation  by  steering  nulls  onto  the  coherent  sources,  but  rather  by  the  weight  phasing  itself.  The 
array  output  was  driven  down  to  receiver  noise  level  in  all  three  cases.  The  plots  for  90°  phase  are  very 
similar  to  what  one  would  ohtain  if  all  three  sources  were  non-coherent,  i.e.,  cancellation  is  achieved  by 
adaptive  null  steering,  in  this  instance. 

Such  severe  sensitivity  to  coherent  source  phasing  in  the  mainbeam  region  produces  different  dis¬ 
tortions  in  tracking  estimates  from  adaptive  I  (sum)  and  A  (difference)  patterns,  as  shown  in  Fig.  14. 
Appendix  C  contains  the  equation  development  for  this  type  of  plot;  however,  the  main  point  here  is  to 
show  the  considerable  changes  in  track  angle  estimates  due  to  phase  variation.  Once  again,  if  all  three 
sources  were  non-coherent,  the  distortion  plot  would  be  stable  and  very  similar  to  the  one  shown  for 
the  90°  phase. 
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Hg  14  —  I  rack  estimate  distortion  resulting  Irom  adaptive  >.  and 
A  tracking  beams  for  coherent  interference  in  the  mainbeam 
region,  same  case  as  Fig.  13 


All-Digital  Tracking  System  Concept 

The  separate  estimation  of  interference  source  data  (total  number,  power  levels,  location  angles, 
coherence)  and  its  utilization  for  improving  the  output  SNR  of  desired  signai  detections  is  a  mode  of 
system  operation  that  has  been  addressed  in  the  literature  a  number  of  times  for  various  applications 
[5,6,15,16].  In  this  section,  we  briefly  review  such  a  system  wherein  the  estimated  data  is  used  to  drive 
a  fully  adaptive  tracking  processor  [27],  The  concept  is  illustrated  in  Fig.  15.  Starting  on  the  left-hand 
side,  the  system  continuously  computes/updates  a  sample  covariance  matrix  R.  Of  particular  signifi¬ 
cance  is  that  R  may  be  dimensioned  either  equal  to  or  less  than  .he  total  number  of  array  elements, 
i.e.,  the  model  order  of  the  estimate  is  selectable  per  subaperture  averaging  option  choice.  Off-line  pro¬ 
cessing  on  R  is  then  conducted  at  periodic  intervals  to  estimate  the  locations  and  relative  power  levels 
of  interference  sources  via  the  most  appropriate  spectral  estimation  algorithms.  The  central  processor 
unit  (CPU)  then  applies  these  data  to  the  computation  of  optimized  adaptive  spatial  filter  weights  for 
the  right-hand  side  of  Fig.  15.  Separation  of  source  estimation  from  adaptive  filter  weight  computation 
can  be  done  accurately  only  in  an  all-digital  processing  system,  but  it  permits  the  following  benefits: 

a.  Estimation  of  coherent  interference  source  locations  for  deliberate  adaptive  null  filter  place¬ 
ment. 

b.  Remembering  slow-changing  or  time-gated  sources,  and  colored-noise  distributions. 

c.  Anticipating  sources  from,  apriori  data  inputs. 
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Fig.  1 3  —  All  digil.il  adaptive  array  tracking  system  concept 


d.  Flexibility  in  time-domain  control  of  the  filtering  to  counter  interference  time  strategies. 

e.  Tracking/cataloging/ranking  sources. 

f.  Efficient  assignment  of  available  DOF. 

g.  Compatible  with  fast-response  adaptive  algorithms,  i.e.,  parallel  algorithm  processing. 

The  right-hand  side  of  Fig.  15  indicates  a  fast-memory  storage  capability  that  is  intended  to  permit 
selected  time  delays  of  the  snapshots  for  feeding  into  the  filter  weights.  The  idea  is  to  synchronize 
selected  snapshots  with  their  filter  weight  updates,  if  possible. 

Finally,  the  filtered  signal  output  residue  is  fed  into  a  bcamformcr  which  is  weighted  to  produce 
the  desired  search  and  monopulse  track  beams  for  target  detection  and  tracking.  The  algorithms  of 
Davis  et  at.,  122],  may  be  applied  for  estimating  the  target  signal  angle  of  arrival,  based  upon  the  out¬ 
puts  of  adaptively  distorted  sum  and  difference  beams.  Appendix  C  discusses  the  equivalence  of  such 
beams  to  the  concept  shown  in  Fig.  15. 

As  an  example,  let  us  apply  this  concept  to  the  coherent  source  case  used  in  Figs.  13  and  14 
wherein  we  would  use  a  16-clement  linear  array  feeding  into  our  all-digital  processor.  An  appropriate 
estimation  algorithm  is  that  of  forward-backward  subaperture  spatial  smoothing  [5,28,29)  combined 
with  eigenanalysis  Appendix  D  describes  the  rudiments  of  this  algorithm,  and  the  results  are  plotted 
in  Fig.  16  in  comparison  with  a  scanned  mainbeam  output.  From  th  s  source  estimation  data,  we  can 
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Fig.  16  —  Comparison  of  niainbcam  scan  vs  spatial  smoothing  processing  for  coherent  source  case 
of  Fig  13,  PEUS  eigenanalysis,  156  snapshots  per  trial 


construct  an  equivalent  covariance  matrix  dimensioned  for  the  full  aperture  per  the  procedure  given  in 
Appendix  A,  and  wc  can  compute  its  inverse  for  obtaining  the  adaptive  filtering.  If  we  define  the  con¬ 
structed  covariance  matrix  as  M,  then  its  inverse  may  be  viewed  as  a  matrix  set  of  adaptive  beantformer 
filter  weights  to  give  us  the  filtered  output  nth  snapshot  vector  E /(n), 

E>(n)  -  E'(n)  (11) 

Conventional  beam  weighting  S*  can  then  be  applied  to  the  filtered  output  residue  to  obtain  the  finai 
output  for  the  nth  snapshot, 

y0(»)-E;(n)S*-E,(n)M'lS*  (12) 

or 

r0(n)  -  E'(n)W0 

where  W0  is  the  familiar  optimum  Wiener  filter  weight. 

Note  that  the  constructed  covariance  matrix,  M,  permits  options  such  as  adding  synthetic  sources 
or  changing  power  levels.  Furthermore,  since  it  is  always  Toeplitt,  solutions  may  be  simplified  some¬ 
what. 


For  the  current  example,  the  computed  adaptive  characteristics  would  be  very  similar  to  those 
plotted  in  Figs.  13  and  14  for  the  90°  phase  angle.  Other  examples  along  with  a  more  detailed  discus¬ 
sion  of  the  processing  may  be  found  in  (27). 

CONCLUSIONS 

Two  conceptual  application  areas  have  been  presented  for  using  spectral  estimation  techniques; 
partially  adaptive  lew-sidelobe  arrays,  and  fully-adaptive  tracking  arrays.  In  both  cases,  improved  spec¬ 
tral  estimation  techniques  are  used  separately  to  acquire  information  about  the  interference  enviren- 
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mcnt  which  is  beyond  that  ordinarily  available  in  a  conventional  adaptive  array.  Examples  discussed 
included  supcrresalution  effects,  relative  power  level  determination,  estimation  of  coherent  sources,  and 
the  tracking/cataloging/ ranking  of  sources.  For  the  partially-adaplivc  area,  the  information  was  used 
for  efficient  assignment  of  a  limited  number  of  DOF  in  a  beamspacc  constrained  adaptive  system  to 
obtain  the  following  benefits  (as  compared  to  a  fully  adaptive  array):  retention  of  low  sideiobes  plus  a 
stable  mainbeam;  considerably  faster  adaptive  response;  reduction  in  overall  cost,  and  greater  flexibility. 
On  the  negative  side,  we  incur  the  risk  of  possible  inferior  cancellation  performance  if  the  interference 
source  situation  is  not  adequately  covered  by  the  assigned  DOF. 

For  the  fully  adaptive  tracking  array  area,  the  information  is  used  in  an  all-digital  processing  sys¬ 
tem  to  obtain  the  benefits  of  stable  nulling  of  coherent  interference  sources  in  the  mainbeam  region, 
efficient  assignment  of  the  available  DOF,  and  a  far  greater  flexibility  in  the  time-domain  control  of 
adaptive  filtering  strategy. 
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Appendix  A 

SNAPSHOT  SIGNAL  MODEL 

Consider  a  simple  linear  array  of  K  elements  as  shown  in  Fig.  Al.  The  received  signal  samples 
are  correlated  in  both  space  and  time,  giving  rise  to  a  two-dimensional  data  problem,  but  we  convert 
this  to  spatial  domain  only  by  assuming  that  narrow-band  filtering  precedes  our  spatial  domain  process¬ 
ing.  Bandwidth  can  be  handled  when  necessary  via  a  spectral  line  approach  [Al]  or  tapped  delay  lines 
at  each  element  [A2],  but  we  did  not  consider  such  extra  compf  ation  essential  to  the  basic  purposes  of 
this  analysis.  Thus,  the  postulated  signal  environment  on  any  given  observation  consists  of  1  narrow- 
band  plane  waves  arriving  from  distinct  directions  9,.  The  RF  phase  at  the  Arth  antenna  element  as  a 
result  of  the  ith  source  would  be  the  product  w,  Xk,  where  Xk  is  the  location  of  the  element  phase 
center  with  respect  to  the  midpoint  of  the  array  in  wavelengths,  and  <u,  is  defined  as 

u>i  —  2ir  sin  9, .  (Al) 


Fig.  Ai  —  Geometry  of  linear  array  and  signal  wavefront 


This  notation  is  deliberately  chosen  to  have  the  spatial  domain  dual  of  sampling  in  the  time  domain,  so 
that  the  reader  may  readily  relate  to  the  more  familial  spectral  analysis  variables.  Sin  9,  is  the  dual  of  a 
sinusoid  frequency  /,,  and  the  Xk  locations  are  the  dual  of  time  sampling  instants  tk.  Note  that  if  our 
elements  are  equally  spaced  by  a  distance  d,  then  Xk  may  be  written, 


d 

k  *+,r 

X 

i  2 

(A2) 


where  X  is  the  common  RF  wavelength.  The  ratio  d I  X  becomes  the  dual  of  the  sampling  time  T  with 
the  cutoff  frequency  equal  to  half  of  the  reciprocal. 
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The  complex  amplitude  of  the  rth  source  at  the  array  midpoint  phase  center  is  p, ,  such  that  we  can 
now  express  the  nth  time-sampled  signal  at  the  Arth  element  as, 

/ 

Ek(n)  —  r)k(n)  +  p,(n)gK(0,)  exp  (j<a,Xk)  (A3) 

/-i 

where  (0,)  is  the  element  pattern  response  in  the  direction  B ,,  and  tjaT/j)  is  the  «th  sample  from  the 
Ar th  element  independent  Gaussian  receiver  noise.  (The  receiver  noise  component  is  assumed  to  be  a 
random  process  with  respect  to  both  the  time  index  n  and  the  element  index  k .)  Equation  (A3)  per¬ 
mits  us  to  construct  a  convenient  column  vector  of  observed  data  in  the  form, 

E(rt)  =  Vp(n)  +  ■»)(«)  (A4) 
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where  V  is  a  K  x  /  matrix  containing  a  column  vector  v,  for  each  of  the  /  source  directions;  i.e., 

vlk  “  exp  (jco,Xk) .  (A5) 

Note  that  Eq.  (A4)  separates  out  the  basic  variables  of  source  direction  in  the  direction  matrix  V, 
source  baseband  signal  in  the  column  vector  p(/i),  and  element  receiver  channel  noise  in  the  column 
vector  Tj(n).  The  vector  E(.»t)  is  defined  as  the  nth  snapshot,  i.e.,  a  simultaneous  signal  sampling 
across  all  A" -array  elements  at  the  nth  time  instant.  These  snapshots  would  nominally  occur  at  the 
Nyquist  sampling  rate  corresponding  to  our  receiver  bandwidth  [A3],  so  that  a  radar-oriented  person 
may  view  them  as  range  bin  time  samplings.  However,  for  source  estimation  purposes,  they  need  not 
necessarily  be  chosen  from  contiguous  range  bins,  and  for  most  applications  it  would  be  highly  desir¬ 
able  to  selectively  lime  gate  the  snapshots  used  for  source  estimation.  For  this  simple  analysis,  let  us 
postulate  that  the  snapshots  are  gated  at  more  or  less  arbitrary  instants  of  time. 

Over  typical  processing  intervals,  the  directions  of  arrival  will  not  change  significantly,  so  that  V  is 
a  slowly  changing  matrix.  In  contrast,  the  signals  p,(n )  will  generally  vary  rapidly  with  time,  often 
unprediclably,  such  that  wc  must  work  with  their  statistical  descriptions.  It  is  assumed  that  the  signals 
are  uncorrelated  with  receiver  noise.  Proceeding  then  from  Eq.  (A4),  we  can  obtiin  the  covariance 
matrix  R  via  application  of  the  expected  value  operator,  #,  or  ensemble  average, 

R  -  £[E  («)£*'(«)]  (A6) 

R  =  VPVv  +  N  (A7) 

where  N  —  &'(ij(n)ij(«) "'],  P  -  &’(p  in) p  in)’1 1,  *  is  the  complex  conjugate,  and  t  is  the  transpose.  N 
is  a  simple  diagonal  matrix  consisting  of  the  receiver  channel  noise  power  levels.  The  diagonal  ele¬ 
ments  of  P  represent  the  ensemble  average  power  levels  of  the  various  signal  sources,  and  off-diagonal 
elements  can  be  nonzero  if  any  correlation  exists  between  the  sources.  Note  that  correlated  far-field 
signals  can  easily  arise  if  significant  specuiar  reflection  or  diffraction  multipavh  is  present. 
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Appendix  B 

EIGENVALUE/EIGENVECTOR  DECOMPOSITION 


When  a  signal  is  known  to  consist  of  pure  sinusoids  in  white  noise,  an  appropriate  procedure  to 
find  the  unknown  frequencies  and  powers  is  the  Pisarenko  spectral-decomposition  procedure  [B 1  ] . 
Although  Pisarenko’s  method  per  se  has  not  found  widespread  use,  it  has  provided  a  fundamental 
eigenanalysis  basis  for  several  closely  related  techniques  which  have  demonstrated  excellent  perfor¬ 
mance.  Among  these  are  the  algorithms  described  by  Reddi  [B2],  the  MUSIC  algorithm  of  Schmidt 
[B3],  the  work  of  Bienvenu  and  Kcpp  [&4],  the  singular  value  decomposition  (or  principal  eigenvector) 
methods  of  Kumaresan  and  Tufts  [B5,B6],  the  eigenassisted  method  of  Evan  et  al.  [B7],  and  the  alge¬ 
braic  approach  of  Bronez  and  Cadzow  [B8]. 

A  key  principle  in  all  of  these  techniques  is  the  geometric  vector  space  relationships  between  the 
spatial  source  vectors  and  the  eigenvectors  of  the  sample  covariance  matrix;  so  we  begin  our  discussion 
on  this  point.  From  the  theory  of  matrices,  we  know  that  a  positive  definite  Hermitian  matrix  such  as 
R  of  Eq.  (A7),  can  be  diagonalized  by  a  nonsingular  orthonormal  modal  matrix  transformation  which 
shall  be  defined  as  the  matrix  Q.  Furthermore,  we  know  that  the  resulting  diagonal  components  are 
the  eigenvalues  of  matrix  R.  In  accordance  with  the  usual  eigenvalue  problem  statements, 

|R  —  /3,?I  I  =  0  anu  Re,  =  /3,2e, ,  (Bl) 

the  /3,?  are  the  eigenvalues  (real  positive  numbers)  of  R,  I  is  the  identity  matrix,  and  e,  are  the  associ¬ 
ated  eigenvectors  These  eigenvectors,  which  are  normalized  to  unit  Hermitian  length  and  are  orthog¬ 
onal  to  one  another,  make  up  the  columns  of  the  Q  matrix. 
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Diagonalizalion  of  R  by  the  Q  matrix  transformation  per  Eq.  (Bl)  may  be  written, 


8 1  0  0 
0  @2  0 


(Q’'RQl  «  1/3, J5,,] 


0  0  8i  ■  ■ 


(B3) 


where  S y  is  the  Kronecker  delta  symbol.  One  can  readily  show  a  construction  of  R  from  its  orthonor¬ 
mal  components, 

*  «  Qfe/6  JQ*' -  £/3A2e*e;'.  (B4) 

A- 1 


Next,  we  introduce  the  term  “principal  eigenvector"  (PE)  to  mean  those  eigenvectors  which 
correspond  to  the  unique  eigenvalues  generated  by  the  spatial  source  distribution;  and  the  term  "noise 
eigenvector"  to  mean  those  eigenvectors  which  correspond  to  the  small  noise  eigenvalues  generated  by 
the  receiver  channel  noise  in  Eq.  (A7).  Under  ideal  conditions,  the  noise  eigenvalues  are  all  identical 
and  equal  to  receiver  channel  noise  power  level  8$ ,  such  that  we  can  factor  Eq.  (B4)  to  emphasize  the 
PE, 

R-  0O1  (B5) 

i-i 

where  q  is  the  number  of  PE.  Comparing  Eq.  (B5)  with  Eq.  (A7)  we  note  that  the  noise  diagonal 
matrices  are  equal,  i.e., 

8$l  -  N  (B6) 

so  that  one  may  equate  the  source  direction  vectors  with  the  PE, 

VPV*'  -  £  I  p,  IV/'  -  £  (8}  -  /f02)c,e/  (B7) 

i-i  i-i 

where  the  ip,  I2  represents  the  expected  power  levels  of  uncorrelated  sources.  Equation  (B7)  embodies 
the  key  principle  that  the  PE  are  linear  combinations  of  the  source  direction  vectors  and  vice-versa.  In 
geometrical  language,  the  v,  define  an  /  dimensional  source  vector  space,  and  the  principal  e,  span  that 
same  vector  space.  Furthermore,  since  the  noise  eigenvectors  are  always  orthogonal  to  the  PE,  then  it 
follows  that  the  noise  eigenvectors  must  occupy  a  subspace  which  is  orthogonal  to  the  source  vector 
space.  To  put  it  another  way,  if  the  noise  eigenvectors  are  viewed  as  antenna  array  element  weights, 
then  they  should  have  pattern  nulls  at  source  direction  angles  because  of  their  orthogonality.  (Note  the 
vivid  demonstration  of  this  point  in  Fig.  3(b).)  Despite  the  fact  that  Eq.  (B7)  is  based  on  ideal  assump¬ 
tions,  it  turns  out  to  be  a  valuable  concept  for  formulating  algorithms,  perhaps  because  it  is  inherently  a 
noise-subtracted  relationship,  and  the  estimates  of  the  PE  are  rather  robust. 


27 


W.  F.  GABRIEL 


When  working  with  finite  sets  of  data  snapshots  which  are  not  ideal,  a  nontrivial  problem  area 
arises  in  determining  which  eigenvectors  to  designate  as  principal,  and  which  ones  result  from  noise. 
This  important  problem  will  be  addressed  after  discussing  the  associated  algorithms. 

Eigcnanalysis  of  Three  Algorithms 

Eigenvalue/eigenvector  decomposition  is  now  applied  to  three  different  spatial  spectrum  estima¬ 
tion  algorithms: 

•  the  MLM  (Maximum  Likelihood  Method), 

•  the  MUSIC  (Multiple  Signal  Classification),  and 

•  the  PEGS  (Principal  Eigenvector  Gram-Schmidt). 

Since  these  algorithms  are  dealt  with  in  a  very  abbreviated  manner  in  the  report,  the  reader  is 
encouraged  to  consult  the  references  given  for  a  better  description  and  understanding  of  the  techniques 
involved. 

The  MLM  (Maximum  Likelihood  Method) 


The  maximum  likelihood  spectral  estimate  is  defined  as  a  filter  designed  to  pass  the  power  in  a 
narrow-band  about  the  signal  frequency  of  interest  and  to  minimize  or  reject  all  other  frequency  com¬ 
ponents  in  an  optimal  manner  [B9,B10],  This  is  identical  to  the  use  of  a  zero-order  mainbeam  direc¬ 
tional  gain  constraint  in  adaptive  arrays  [Bill,  where  the  spatial  spectrum  would  be  estimated  by  the 
output  residua)  power  PG  from  the  optimized  adapted  array  weights, 

P0(f»  =  W0*'RW0.  (B8) 


Where  W0  is  the  optimum  adaptive  Wiener  filter  weight,  and 

W0«MR~!S*  (B9) 


where  S*  is  the  usual  mainbeam  weight  vector  for  steering  angle  9 ,  and  /u  is  a  complex  number. 
Under  the  zero-order  gain  constraint,  we  require  S'W0  *=  1,  whereupon  m  becomes 

,i  -  [S'R-’S*!"1 .  (BIO) 

Substituting  m  and  W0  into  Eq.  (B8)  results  in 


P0(9) 


I 

S'R-'S* 


(BID 


Upon  sweeping  the  steering  vector  S*  for  a  given  covariance  matrix  inverse,  P^iO)  estimates  the  spatial 
spectrum. 


In  terms  of  eigenvalue/eigenvector  decomposition,  we  can  take  the  inverse  of  Eq.  (B3)  and 
express  R_1  in  the  form. 
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(B12) 


Here,  we  see  that  this  older  algorithm  simply  uses  all  of  the  eigenvalues/eigenvectors.  One  advantage 
in  this  decomposition  is  that  Eq.  (B12)  can  be  substituted  into  Eq.  (Bll)  to  form  a  simple  summation 
of  eigenvector  beams  referenced  to  the  receiver  noise  power  level. 
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where  gk(0)  =  S'e*.  This  permits  an  insight  into  the  peak  values  which  occur  at  the  nulls  of  the  noise 
eigenvector  beams;  i.e.,  we  get  an  evaluation  of  relative  source  power  level  if  the  source  location  is 
resolved. 

The  MUSIC  (Multiple  Signal  Classification) 

This  algorithm  was  suggested  by  Schmidt  [B3]  to  provide  asymptotically  unbiased  estimates  of  the 
number  of  signal  sources,  directions  of  arrival,  strengths  and  cross  correlations  among  the  directional 
waveforms,  polarizations,  and  strength  of  noise/interference.  His  geometrical  vector  space  description 
and  interpretation  is  clearly  presented  and  was  used  as  the  basis  for  the  one  above.  Essentially,  this 
MUSIC  algorithm  selects  and  uses  only  the  noise  eigenvectors  to  solve  for  the  directions  of  arrival. 
This  is  tantamount  to  approximating  R_1  in  Eq.  (BID  by  the  noise  eigenvectors  only;  i.e., 

let  R“'  ~  X  TT  e<e<‘'  (B14) 

<-«+i  \  j 

where  q  is  the  number  of  principal  eigenvectors.  The  same  indexing  would  apply  as  noise  eigenvector 
beams  in  Eq.  (B13),  where  we  note  that  the  ratio  of  eigenvalues  would  now  become  unity  (or  close  to 
it). 


This  algorithm  does  indeed  produce  very  large  peaks  in  Pn(0)  for  good  covariance  matrix  esti¬ 
mates,  because  of  the  aforementioned  orthogonality  of  the  noise  eigenvectors  to  the  source  vector 
space.  Its  performance  is  usually  far  superior  to  the  older  MLM  algorithm  in  resolving  closely  spaced 
source  directions.  In  addition,  Schmidt  points  out  that  once  the  directions  of  arrival  have  been  found, 
the  direction  matrix  V,  in  Eq.  (A4)  and  (A7)  becomes  available  and  may  be  used  to  compute  the 
sou  ice  power  matrix  F.  We  form  the  special  matrix  *J, 

U  -  IV'VJ-'V'  (B15) 

and 


UVPV'lT'  -  P. 

If  the  matrix  Li  exists  or  can  be  closely  approximated,  then  we  can  app'y  it  to  the  noise-subtracted 
covariance  matrix  of  Eq.  (A7)  to  solve  for  P;  i.e., 

P  “  U(R  —  N)ll"'  (B16) 


P  = 


P\*P\  P\  *Pi  ...  P]'Pi 


P2*P  1  Pl*P2  ■  •  P2*P, 


PimP\  Pi  'Pi  P;  *Pi 


Note  that  the  diagonal  elements  of  P  represent  power  estimates  of  the  sources,  and  that  the  nonzero 
off-diagonal  elements  represent  estimates  of  the  correlations  existing  between  partially  coherent 
sources. 
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The  ability  to  solve  for  the  power  estimates  is  of  great  importance  in  distinguishing  "false  alarms" 
and  in  selecting  the  sources  of  interest.  We  recommend  References  B7,  B6,  and  B8  which  at  :  either 
related  to  or  give  a  comparative  analysis  of  the  MUSIC  algorithm. 

The  PEGS  (Principal  Eigenvector  Gram-Schmidt) 

Several  eigenvalue/eigenvector  decomposition  techniques  described  in  References  B7,  B2,  B6,  ar.d 
B8  are  based  on  the  principal  eigenvectors  (PE)  with  some  type  of  constraint  imposed  on  the  optimum 
weight  vector. 

This  subgroup  of  PE  methods  is  of  interest  in  the  current  work  because  of  their  generally  superior 
performance  characteristics.  An  intuitive  reasoning  behind  their  use  is  that  the  estimates  of  the  PE  are 
robust;  i.e.,  they  tend  to  remain  relatively  stable  from  one  data  record  to  the  next,  whereas  the  noise 
eigenvectors  tend  to  fluctuate  because  of  noise  perturbations.  In  addition,  the  PE  methods  are 
inherently  a  noise  subtraction  technique  similar  to  noise  power  cancellation  algorithms  which  attempt  to 
remove  the  noise  bias  term  that  appears  along  the  main  diagonal  of  the  covariance  matrix. 


Let  us  begin  development  of  our  PEGS  algorithm  by  decomposing  the  inverse  of  the  covariance 
matrix  as  given  in  Eq.  (B12),  normalized  by  receiver  noise  power, 


\  &  ~  flo 

l  fi?  J 


o  o 


(B 1 7) 


where: 

«,  -  e/'S* , 
p.'  -  p. / /3q ,  and 

S*  «=  a  quiescent  array  weight  vector. 


In  the  limit  of  noise-free  data,  Eq.  (B18)  is  suggestive  of  a  simple  Gram-Schmidt  vector  subtrac¬ 
tion  from  S*  in  which  we  would  form  an  optimum  weight  that  would  be  orthogonal  to  the  PE,  and 
therefore,  orthogonal  to  the  source  direction  vectors  per  Eq.  (B7).  Thus,  let  us  formulate  a  PEGS  algo¬ 
rithm  by  defining  the  optimum  weight  Wf  from  Eq.  (B 1 8)  as 

W,- S' (B19) 

i-i 

where: 

a.-e.^SV 

W,  possesses  the  necessary  orthogonality  to  noise-free  source  direction  vectors, 

<  Wf ,  v,  >  =  0  (B20) 
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and  may  readily  incorporate  the  option  of  unit  Hermitian  length  if  desired, 

|WPP=  1.  (B21 ) 

The  PEGS  algorithm  as  applied  in  this  report  used  an  end  clement  weighting  for  S*;  i.c., 

S'-  [0  0,  ...,0.11'.  (B22) 

Culling  Principal  Eigenvalues 

The  number  of  principal  eigenvalues  is  usually  directly  related  to  the  number  of  sources  which,  in 
practice,  are  not  known  and  must  be  estimated.  One  of  the  early  estimation  techniques  which  has  often 
been  used  is  the  A1C  (Akaike  Information  Criterion)  [B 1 2.B  1 3)  This  criterion  has  been  successfully 
applied  to  many  model  identification  problems  in  engineering  and  statistics,  including  the  well  known 
problem  of  determining  the  order  of  an  autoregressive  (AR)  process  [B 13).  Recent  work  reported  by 
Wax  and  Kailath  (B!4)  presents  a  new  approach  based  on  the  A1C,  which  eliminates  the  need  for  any 
subjective  judgment  in  the  decision  process;  i.e.,  the  procedure  does  not  require  any  subjectively  chosen 
threshold.  This  new  approach  was  implemented  during  the  current  investigation  and  was  found  to  be 
very  effective  for  most  of  the  examples  tested.  In  addition  to  the  Wax-Kailath  A1C  approach,  we  also 
used  a  second  effective  technique  which  is  based  on  the  following  three  processing  operations: 

•  An  initial  decreasing-magnitude  sort, 

«  Culling  per  coarse  magnitude  threshold,  and 

•  Culling  per  sensitive  threshold  based  on  a  quadratic  curvefit  predictor. 

The  entire  procedure  is  listed  in  Ref.  (B15]  as  a  Fortran  IV  computer  code. 
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Appendix  C 

TRACKING  BEAMS  AND  ADAPTIVE  DISTORTION 


The  tracking  beams  used  in  Section,  An  Adaptive  Art  ay  Tracking  Application,  are  based  on  the 
selection  of  an  adjacent  pair  of  orthogonal  uniform  illumination  beams  generated  by  a  Butler  matrix 
bcamformer  [Cl]  transformation.  The  transformation  matrix  B  for  a  linear  array  with  half-wavelength 
element  spacing  will  have  individual  matrix  components  of  the  form, 
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K  +  1 
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(Cl) 


where: 


m  =  beam  index 


k  =  element  index 
K  =  total  number  of  elements. 
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The  beamformer  output  vector  E  is  expressed, 

E  -  B'E  (C3) 

Figure  Cl  (a)  illustrates  typical  sinx/x  patterns  of  an  adjacent  pair  of  beams  for  our  16-element  array 
example.  The  familiar  sum  (£)  and  difference  (A)  tracking  beam  outputs  are  then  obtained  from  this 
adjacent  pair  via  the  3  oB  hybrid  junction  as  shown  in  Fig.  C2,  where  we  note  that  the  A  beam  is  in 
quadrature  phase  relationship  to  the  I  beam.  Expressed  in  terms  of  element  weights,  the  I  beam 
weight  vector  S  and  the  A  beam  weight  vector  D  may  be  written  via  Eq.  (C2), 

S-  j  lbm+)  +  bml 

D-jlb^,  — bj.  (C4) 

The  uniform  illumination  vectors  bm  and  bm+i  result  in  ccsine  illuminations  for  S  and  D,  which  are 
shown  plotted  in  Fig.  Cl(b)  for  our  16-element  linear  array  example. 

Monopulse  tracking  [C2]  involves  an  angle  estimate  for  each  pulse  (snapshot)  containing  the  tar¬ 
get,  and  it  is  computed  from  the  ratio  of  A/E.  From  the  cosine  illumination  beams  of  Fig.  Cl,  note 
that  wc  can  form  the  approximation. 
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8  is  the  angle  of  the  target  from  borcsight, 

B  is  the  angle  to  the  first  null  of  the  sum  pattern,  and 
C  is  a  constant  which  depends  upon  the  particular  illuminations. 


(C5) 


BEAMFORMEN 

OUTPUTS 


I  A 

BEAM  BEAM 

Fig.  C2  —  Block  diagram  of  (racking  beam  formation 
via  Butler  matrix  bcamformer 

For  our  example,  C  *=  3.8"  and  B  ~  10.8°.  Thus,  given  values  of  A/Z,  we  can  solve  for  track  angle 
estimates,  8/fl,  from  Eq.  (C5). 

Next,  let  us  address  track  estimate  distortion.  Whenever  one  performs  spatial  filtering  as 
described  in  Section,  All-Digital  Tracking  System  Concept,  a  distortion  of  received  plane  wavefronts 
occuts  because  the  spatial  insertion  loss  generally  vanes  as  a  function  of  sin  9. 

This  problem  was  first  addressed  in  the  literature  by  Davis,  Brennan,  and  Reed  lC3)  who  pro¬ 
posed  an  algorithm  for  estimating  the  angle  of  arrival,  based  on  the  outputs  of  adaptively  distorted  sum 
and  difference  beams.  They  used  approximations  to  the  optimum  angle  estimator  which  permitted 
correction  of  distortion  at  the  tracking  beam  boresight  position,  and  they  demonstrated  good  perfor¬ 
mance  via  simulation  for  sidclobe  and/or  mainbeam  interference. 

In  an  all-digital  system,  we  know  our  adaptive  filter  weights.  Therefore,  we  can  compute  the 
resultant  distortion  error  throughout  the  tracking  beam  region.  In  Section,  All-Digital  Tracking  System 
Concept,  above,  wc  showed  that  it  makes  no  difference  whether  we  apply  our  quiescent  beam  weights 
to  the  spatially  filtered  signals,  or  the  equivalent  adaptive  weights  to  the  unfiltered  signals.  Thus,  we 
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may  apply  the  monopuise  sum  and  difference  weights  of  Eq.  (C4)  to  the  spatially  filtered  output  resi¬ 
due  signal  vector  of  Eqs.  (11,  12)  in  Section,  All-Digital  Tracking  System  Concept,  and  obtain  the 
equivalent  forms. 


I  -  S'E,  -  fi'M  ‘S  -  E'WS 
and 

A  -  D'E,  -  E'M  ’D-  E'Wa  , 


(C6) 

(Cl) 


where  Ws  and  XV  d  are  now  the  equivalent  adapted  (and  distorted)  sum  and  difference  beam  weights. 
The  distorted  ratio  A/E  can  be  computed  for  any  d  ection  vector  E,  thus  giving  us  the  distortion  curve 
across  the  entire  tracking  angle  region. 
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Appendix  D 

FORWARD-BACKWARD  SUBAPERTURE  AVERAGING 


This  is  an  excellent  technique  for  increasing  the  effective  averaging  of  our  sample  covariance 
matrix  SCM  when  needed,  and  it  may  readily  be  implemented  if  the  antenna  array  elements  are  identi¬ 
cal  and  equ.-.lly  spaced.  Figure  D1  illustrates  this  technique.  We  form  a  reduced  dimension  subaperture 
of  L.  elements,  where  L  must  be  less  than  ihe  total  number  of  array  elements  K.  Starting  from  the 
left-hand  siae,  the  subaperture  samples  its  first  snapshot  as  elements  1  through  I,  then  bumps  to  the 
right  by  1  and  samples  its  second  snapshot  as  elements  2  through  ( L  +  1),  then  bumps  to  the  right  by 
1  and  samples  its  third  snapshot  as  elements  3  through  ( L  +  2),  etc.  After  bumping  across  to  the  ATh 
element,  we  will  have  accumulated  {K  —  L  1)  subaperture  snapshots  from  one  overall  array  data 
snapshot,  such  that  we  can  increase  our  SCM  averaging  by  that  same  factor.  This  subaperture  motion 
from  left  to  right  produces  what  is  generally  termed  "forward  averaging."  The  technique  may  be  applied 
to  any  SCM  method.  For  example,  the  simple  block  averaging  of  Eq.  (2)  in  section,  Low-Sidelobe 
Eigenvector  Constraint  becomes 

R/  -  ~s(k-TT 17  i  Ieu,/)E(«,/)*-i  (dd 

where  E(n.r)'  *=  [E,(n),  El+i(n),  ...  ,  and  R/-  is  the  new  reduced  L  x  L  dimension  SCM. 

Note  that  E  («,/)  may  be  expressed  as  the  matrix  product, 

E  {n ,/)  —  l,E  in)  (D2) 

where  I,  is  a  special  L  x  K  rectangular  sampling  matrix  in  which  the  /  index  denotes  the  first  column 
where  the  identity  matrix  I  begins.  For  example,  the  I,  matrix  for  L  =  3,  K  *=  6,  and  /  =  2  would  be 

I,  -  [O  j  I  j  O]  (D3) 
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Fig.  D1  —  Forward-backward  shift  movement  for  a  reduced 
dimension  sampling  subapenure  where  L.  -  4,  along  a  linear 
antenna  array  of  K  -  15  elements 


An  I,  matrix  may  be  used  to  reduce  the  number  of  computations  by  multiplying  the  SCM  to  obtain 

]  (K-L+U 


W  F  G  A  BRILL 


Equations  (D4)  and  (D!)  give  similar  results,  and  both  are  mathematical  expressions  of  the  additional 
spatial  averaging  or  "smoothing"  that  is  obtained  via  the  moving  subaperture  technique. 

Furthermore,  the  averaging  can  be  doubled  again  by  reversing  our  subaperture  at  the  right-hand 
side  and  bumping  across  to  the  left-hand  side  in  similar  fashion;  however,  it  requires  conjugating  the 
subaj  erture  snapshots.  This  subaperture  motion  from  right  to  left  produces  what  is  generally  termed 
"backward-ave.  aging." 

If  we  define  the  reversed  and  conjugated  vector,  E*,  in  terms  of  our  array  element  data  samples  of 
equation  (A3),  then  the  A th  element  signal  sample  becomes 

£*‘~£V*+1  ;  *=1-2,3 . K  (D5) 

The  reader  can  verify  that  the  resulting  SC'M  will  be  an  index-exchanged  Rj.  and  that  we  can  combine 
the  two  into  a  final  forward-backward  average  SCM  which  is  denoted  as  the  reduced  L  x  L  matrix  R;ft: 

R,-i  =  1/2  [R,  +  R}1 .  (D6) 

Note  that  R/t  is  a  symmetric  matrix,  but  is  generally  not  Toeplilz.  References  D1  and  D2  are  recom¬ 
mended  for  further  detailed  dis  ussion  of  the  technique. 

Although  forward-backward  subaperture  averaging  is  a  very  simple  concept,  it  usually  produces 
remarkable  improvements  in  output  estimates,  and  it  becomes  crucial  to  processing  in  the  following 
situations: 

•  when  oniy  a  few  data  snapshots  are  avaiiabie  per  SCM  computation.  Note  that  the  method  can 
be  used  even  under  the  extreme  condition  of  only  a  single  snapshot. 

•  When  significant  coherence  exists  between  spatial  sources  as  for  example  in  multipath  situa¬ 
tions  involving  a  specular  reflection.  For  this  particular  condition,  the  fields  arriving  at  the 
aperture  are  nonstationary  in  space  and  the  SCM  is  not  Toeplitz  [D1  ,D3,D4,D5] 

A  caveat  concerning  this  averaging  technique  is  that,  as  the  dimension  L  of  the  subaperturr 
becomes  smaller,  the  subaperture  antenna  gain,  resolution,  and  degrees  of  freedom  decrease.  Thus, 
the  advantage  of  increased  averaging  must  always  be  balanced  against  these  factors,  and  it  is  usually 
prudent  to  process  with  as  large  a  subaperture  dimension  as  possible. 
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